Astronomy & Astrophysics manuscript no. 8210Hat 


© ESO 2008 


February 1, 2008 





H-atom bombardment of CO2, HCOOH and CH3CHO containing ices 

S. E. Bisschop, G. W. Fuchs, E. F. van Dishoeck, and H. Linnartz 



Raymond and Beverly Sackler Laboratory for Astrophysics, Leiden Observatory, Leiden University, P.O. Box 9513, 2300 RA Leiden, 
the Netherlands 



o 
o 

(N 
Oh 



43 

Oh 
6 



> 
in 

On 

in 

(N 

ON 

o 
l> 
o 



>< 



Received; accepted 



ABSTRACT 



Context. Hydrogenation reactions are expected to be among the most important surface reactions on interstellar ices. However, solid 
state astrochemical laboratory data on reactions of H-atoms with common interstellar ice constituents are largely lacking. 
Aims. The goal of our laboratory work is to determine whether and how carbon dioxide (C0 2 ), formic acid (HCOOH) and acetalde- 
hyde (CH 3 CHO) react with H-atoms in the solid state at low temperatures and to derive reaction rates and production yields. 
Methods. Pure C0 2 , HCOOH and CH3CHO interstellar ice analogues are bombarded by H-atoms in an ultra-high vacuum experi- 
ment. The experimental conditions are varied systematically. The ices are monitored by reflection absorption infrared spectroscopy 
and the reaction products are detected in the gas phase through temperature programmed desorption. These techniques are used to 
determine the resulting destruction and formation yields as well as the corresponding reaction rates. 

Results. Within the sensitivity of our set-up we conclude that H-atom bombardment of pure CO2 and HCOOH ice does not result in 
detectable reaction products. The upper limits on the reaction rates are <7xl0~ 17 cm 2 s -1 which make it unlikely that these species 
play a major role in the formation of more complex organics in interstellar ices due to reactions with H-atoms. In contrast, CH3CHO 
does react with H-atoms. At most 20% is hydrogenated to ethanol (C2H5OH) and a second reaction route leads to the break-up of 
the C-C bond to form solid state CH4 (~20%) as well as H2CO and CH3OH (15-50%). The methane production yield is expected to 
be equal to the summed yield of H2CO and CH3OH and therefore CH4 most likely evaporates partly after formation due to the high 
exothermicity of the reaction. The reaction rates for CH 3 CHO destruction depend on ice temperature and not on ice thickness. The 
results are discussed in an astrophysical context. 
Conclusions. 

Key words, astrochemistry - molecular data - ISM: molecules - methods: laboratory - molecular processes 
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1. Introduction 

It is generally assumed that at the low temperatures in inter- 
stellar clouds, thermal hydrogenation of molecules on icy grain 
surfaces is the main mechanism to form more complex satu- 
rated species (Tielens & Charnley 1997). This is due to the rel- 
atively high abundance of H-atoms in the interstellar medium 
as well as their high mobility even on cold grains. However, 
most of the reactions in the propos ed reaction sch e mes have no t 
yet been measu r ed exp e rimentally. [Hiraoka et al.l ( 1 1 994 120021 
Watanabe litail (|2004|) , iHidaka et all d2004 and iFuchs etal 
d2007l) have studied reactions of thermal H-atoms with CO ice in 
the laboratory, and shown that H2CO, and at higher fluxes also 
CH3OH, are readily formed at temperatures as low as 12 K. It 
thus seems likely that other species will also be able to react with 
H-atoms to form saturated grain-surface products. These species 
may be the starting point for an even more complex chemistry 
that occurs at higher temperatures or by energetic processing due 
to UV or cosmic rays in the ice. Eventually, the ices will evap- 
orate when heated by a pr otostar, leading to the complex organ - 
ics seen in hot cores (e.g. jBlake etaflfl987llikeda et alll200lh . 
The aim of this paper is to study the reactivity of a number of as- 
trophysically relevant molecules with H-atoms in interstellar ice 
analogues at low temperatures to test the proposed thermal hy- 
drogenation reaction scheme and to characterize which products 
are formed and which mechanism is involved. 
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Interstellar ices contain both simple and complex species 
(see lEhrenfreund et al.L lT999). The most abundant ice molecules 
are H2O, CO and CO2, which have very strong vibrational 
modes. The spectroscopic identification of other less abundant 
ices, such as HCOOH and CH3CHO studied here, relies on 
weaker bands i.e., the OH and CH bending modes, vb(OH/CH), 
of HCOOH at 7.25 /urn and the CH 3 deformation, v D (CH 3 ), of 
CH3CHO at 7.41 yizm. Solid state abundances of HCOOH are 
1-5% in both l ow and high mass s t ar forming regions with 
respect to U 7 Q dSchutte et all fl997l [1991 iGibb et all 12001 
iBoogert et aill2004l) . The detection of CH3C HO is less certain 
but abundances up to 10% have been reported dGibb et all 12004 
iKeaneet all 120011) . 

The specific species studied in this paper are CO2, HCOOH 
and CH3CHO, molecules that take a central place in inter- and 
circumstellar hydrogenation reaction schemes (Fig. [T). CO2, 
HCOOH and CH 2 (OH) 2 differ only in their number of hydro- 
gen atoms. It is therefore possible that they are related through 
successive hydrogenation reactions. Previous laboratory experi- 
ments of H2 and CO2 h ave resulted in th e formation of HCOOH 
on ruthenium surfaces (Qgo et ail 120061) . but this reaction may 
have been mediated by the catalytic surface. 

Another series of organics that are thought to be linked 
through successive hydrogenation are CH2CO, CH3CHO and 
C2H5OH (Fig. [TJ. Ethanol is indeed detected in warm gas 
phase environments in star-forming region s but CH2CO and 
CH3CHO are f ound m ostly in colder gas dlkeda et ail [2001; 
Biss chop et all l2007bl) . This may be either due to very effi- 
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Fig. 1. Potential reaction routes for hydrogenation of CO2 and 
CH2CO ice. Only stable products are shown. 

cient conversion of CH3CHO into C2H5OH or because the lat- 
ter species may be formed through another route. In particular, 
astronomical observations show a constant CH3OH/C2H5OH ra- 
tio which indicates that these two species are chemically linked 
dBisschopetaUl2007bl) . 

H-addition reactions in astrophysically relevant ices have 
been studied previously through UV photolysis experiments, 
where the hydrogen atoms are produced by d issociation of a suit- 
able precursor mo l ecule, o ften H 2 O ice (e.g. , Ewing et al., 1960; 



Millig an& Jacoxl 
Allamandola et al., 



1964 



197 1 ; Van Iizendoorn et al. 



1983 



, .19881: iGerakines et al.L I2000h iMoore et aL 
200lHWuetalll2002h" . Although these experiments give a useful 



indication whether certain hydrogenation reactions may or may 
not occur, their results cannot be compared directly with those 
obtained from the laboratory studies mentioned above, nor can 
they be used to quantitatively test reaction schemes such as those 
in Fig.[T](see § !2.2l for details). 

This paper is organized as follows: §|2] explains the experi- 
mental method, §[3] focuses on the data reduction and analysis, 
§ |4]-[6] discuss the results, derived reaction rates and chemical 
physical mechanisms for hydrogenation reactions with HCOOH, 
CO2 and CH3CHO, §|7]presents the astrophysical implications 
and finally §[8] summarizes the main conclusions of this study. 

2. Experiments 

2.1. Our experiment 

The experiments are performed using a new ultra-high vacuum 
set-up that comprises a main chamber and an atomic line unit. 
The details o f the operati o n and performance of the set-up are 
described by Fuch s et al] (I2007I) . The main chamber contains 
a gold coated copper substrate (2.5x2.5 cm 2 ) that is mounted 
on top of the cold finger of a He cryostat. Temperatures can 
be varied between 12 and 300 K with 0.5 K precision using a 
Lakeshore 340 temperature control unit and are monitored with 
two thermocouples (0.07% Au in Fe versus chromel) that are 
mounted on the substrate face and close to the heater element. 
The typical pressure in the main chamber during operation is 
better than 5xl0~ 10 mbar. 

Pure ices of 13 C 18 2 (97% purity, Icon), 12 C 18 2 (97% pu- 
rity, Icon), HCOOH (98% purity, J. T. Baker) and CH3CHO 
(99% purity, Aldrich) as well as mixed ices of 13 C 18 02 with 
H2O (deionized) and CO (99.997% purity, Praxair) are studied 
(see Table Q] for an overview of the mixture ratios, ice thick- 
nesses and ice temperatures). The two isotopologues of CO2 are 
used to distinguish between atmospheric CO2 and solid CO2 
processed by H-atoms. The ices are grown at 45° with a flow 
of 1.0xl0~ 7 mbar s _1 where 1.3xl0~ 6 mbar s _1 corresponds 
to 1 monolayer (ML) s _1 . The temperatures of the ices range 
from 12 to 20 K and their thicknesses are chosen between 8 and 
60 ML. 

H-a toms are produced in a well-studied thermal-cra cking 
device dTschersich & von Boninl 11998b iTschersi ch. 2000). The 
dissociation rate and resulting H-atom flux depend on the tem- 
perature and pressure which are both kept constant during a 



single experiment. The temperature of the heated Tungsten fil- 
ament, 7V is ~2300 K in all experiments and the H+H2 flow 
through the capillary in the atom line is either l.OxlO -4 or 
1.0xl0~ 5 mbar s _1 . For the latter pressure the calculated dis- 
sociation rate, aa s , in the atomic line is 0.45. The atoms that exit 
the source are hot. Before the H-atoms enter the main chamber, 
they pass through a quartz pipe and equilibrate to ~300 K. The 
minimum number of collisions of H-atoms with the quartz pipe 
is 4 due to the nose-shaped form of the pipe. Since H-atoms are 
thermalized with the surface after only 2-3 collisions, it is ex- 
pected that most H-atoms will have temperatures equal to the 
quartz pipe of 300 K. Due to collisions with the walls of the pipe 
and with each other, a fraction of the atoms recombines to H2. 
The effective dissociation fraction, a??, and H-atom flux on the 
samp le surface are therefore calculated to be lower than given 
bv lTschersichl d2000b . i.e., 0.13 and 5.0xl0 14 cnr 2 s -1 , or 0.20 



and l.OxlO" 5 mbar s -1 , respectively (iFuchs et all 120071) . Note 
also that, the absolute number of H-atoms on the surface at a 
given time is not equal to the surface flux, because processes 
such as scattering and recombination take place on the sur- 
face. Theoretical simulations show that a steady-state H-atom 
coverage of 5.0xl0 14 crrT 2 s with an error of a factor 2 is 
quickly reached in this regime independent of the exact H-atom 
flux (Cuppen, private communication). A more extensive dis- 
cu ssion of the d erivation of these steady-state numbers is given 
by IFuchs et all d2007l) . The time and H-fluence, i.e., the total 
number of atoms crrT 2 integrated over time in each experiment 
are listed in Table Q] At the temperatures and fluxes used in 
our experiment, the substrate surface will be covered with H2 
in a few seconds. Since H2 molecules do not stick to other H2 
molecules, the maximum coverage with H2 will only be a few 
monolayers. H-atoms therefore have to diffuse through the cold 
H2 layer before reaching the ice and will be completely thermal- 
ized with the surface at the moment they encounter the ice sam- 
ple. Experiments with only H2 molecules have been performed 
for comparison by setting the source temperature 7\y to ~600 K 
as an additional check to confirm that reactions are due to H- 
atoms and not to H2 molecules (see Table[T|i. 



The ices are monitored by Reflection Absorption Infrared 
Spectroscopy (RAIRS) using a Fourier Transform Infrared 
Spectrometer, covering 4000-700 cirT 1 with a spectral resolu- 
tion of 4 cm" 1 . The infrared path length is the same for all 
experiments. Typically 512 scans are co-added. An experiment 
starts with a background RAIR scan and subsequently an ice 
is deposited onto the substrate surface. Another RAIR spec- 
trum is taken after deposition to determine the initial number 
of molecules in the ice. An additional background spectrum is 
recorded afterward such that subsequently recorded spectra yield 
difference spectra between ices before and after H-atom bom- 
bardment. The next step in the experiments is the continuous H- 
atom bombardment of the ice during which a RAIR scan is taken 
every 10 minutes. After 3 hrs the experiment is stopped and 
a Temperature Programmed Desorption (TPD) spectrum is ob- 
tained using a quadrupole mass spectrometer. The ramp speed is 
2 K min _I and is continued until the temperature reaches 200 K. 
Control experiments with pure ices of CO2, C2H5OH, CH4 and 
CH3OH are studied to determine RAIR band strengths and to 
calibrate the production yields that are measured using the mass 
spectrometer (see §[33}. In those cases RAIR spectra are taken 
right after deposition and subsequently TPD spectra are recorded 
as described here for the other experiments. 
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Table 1. Summary of all H-atom bombardment experiments for 
pure CO2, HCOOH and CH3CHO ices, as well as the mixed 
morphologies. The total time that the ices are exposed to H- 
atoms are indicated with t and the temperature of the tungsten 
filament with Tw- A value Tw of 2300 K refers to an experiment 
in which H2 is dissociated and of 600 K to control experiments 
of ices bombarded with H2-molecules. 



Thickness 


T 

1 ice 




7\ v 


H-fluence 


(ML) 


(K) 


(min) 


(K) 


(molecules cnT 2 ) 


C0 2 


15 


12.5 


180 


2300 


5.4(18) 


15 


14.5 


180 


2300 


5.4(18) 


15 


14.5 


180 


600 




C0 2 :H 2 


15" 


14.5 


180 


2300 


5.4(18) 


15" 


14.5 


180 


2300 


5.4(18) 


15' 


14.5 


180 


2300 


5.4(18) 


15" 


14.5 


120 


2300 


3.6(18) 


15" 


14.5 


180 


600 




CO: 13 C 18 2 


30" 


14.6 


180 


2300 


5.4(18) 


15'' 


14.6 


180 


2300 


5.4(18) 


30'' 


14.6 


180 


2300 


5.4(18) 


45'' 


14.6 


180 


2300 


5.4(18) 


15'' 


14.6 


180 


2300 


5.4(18) 


30 e 


14.6 


180 


2300 


5.4(18) 


HCOOH 


20 


12.5 


240 


2300 


7.2(18) 


20 


40.0 


240 


2300 


7.2(18) 


20 


12.5 


240 


600 




HCOOHHjO' 


40 


12.5 


180 


2300 


5.4(18) 


CH3CHO 


16.2 


14.5 


180 


2300 


5.4(18) 


7.8 


14.5 


180 


2300 


5.4(18) 


11.4 


14.5 


180 


2300 


5.4(18) 


13.5 


14.5 


180 


2300 


5.4(18) 


18.8 


14.5 


180 


2300 


5.4(18) 


21.2 


14.5 


180 


2300 


5.4(18) 


22.1 


14.5 


180 


2300 


5.4(18) 


56.0 


14.5 


180 


2300 


5.4(18) 


45.8 


14.5 


180 


2300 


5.4(18) 


11.4 


12.4 


180 


2300 


5.4(18) 


11.3 


17.4 


180 


2300 


5.4(18) 


11.2 


19.3 


180 


2300 


5.4(18) 


11.7 


14.5 


180 


600 





"39:61% 13 C 18 0,:H,0, 6 22:78% 13 C 18 2 :H 2 0, c 48:52% 13 C' 8 2 :H 2 0, 
rf 45:55% 13 C 18 2 :CO, e 80:20% 13 C 18 2 :CO, 'HCOOH:H 2 20:80%. 



2.2. Comparison with other hydrogenation experiments 

A large number of photolysis experiments of astrophysically rel- 
evant ices exist where H-atoms are produce d through photo dis- 
sociation of H 2 Q o r othe r precu r sors (e.g.. lEwing et all 1960; 



Milligan & Jacox, 
Allamandola et al. 



1964 



1988 



1983 



1971; Van Ii zendoorn et al 
iGerakines et all 12000b iMoore et all 
2001; lWu et alll2002l) . These give useful information on poten- 
tial hydrogenation reactions schemes, but do not give specific 
information about reaction rates (see for example § 14.31 & l6~3l . 
Also the question whether thermal hydrogenation reactions can 
be responsible for newly formed species is not answered, for sev- 
eral reasons. First, the hydrogen atoms resulting from photolysis 
are produced in situ inside the ice with an excess energy of sev- 



eral eV; such atoms can travel s ignificant distances through the 
ice (e.g., lAndersson et all 120061) and a reaction may take place 
before thermalization is achieved. Thus, activation energy bar- 
riers can be overcome, in contrast to thermal hydrogenation re- 
actions where this is less probable. Second, although dilution in 
an inert Ar matrix can stabilize the H-atoms, the H-atom flux 
on the reactants remains poorly characterized. Third, other re- 
active products such as OH are also formed by photolysis of 
H2O which makes it hard to discriminate the different effects. 
It is not possible to study pure ices through this method since 
H2O or another precursor is always needed to provide a source 
of H-atoms. Finally, the photolysis experiments reported so far 
have been carried out under high vacuum conditions (typically 
10~ 7 mbar) in which several monolayers of background gases 
(mostly H2O) are accreted in less than a minute, providing addi- 
tional molecules that can be photolyzed during the experiments. 
This may affect the outcome. In contrast, our experiments and 
the previously mentioned surface science experiments on CO hy- 
drogenation are performed under ultra-high vacuum conditions 
(typically 10" 10 mbar) in which less than a monolayer of back- 
ground gas (mostly H2) is accreted during the time-scale of the 
experiments (a few hours). Here the H-atoms are formed by a 
microwave source or by thermal cracking and are thermalized to 
room temperature or less before striking the ice surface, rather 
than being produced inside the ice. 



3. Data analysis 

3.1. R AIR analysis 

Different frequency ranges are selected for baseline subtraction 
that depend on the species under study. Fourth order polynomial 
baselines are fitted to the recorded RAIR spectrum. Additionally, 
local third order polynomial baselines are subtracted around the 
features of interest to accurately determine the integrated ab- 
sorption. The frequency ranges are given per ice morphology 
in Table [2] From the integrated intensity of the infrared bands, 
the column density of species X in the ice is calculated through 
a modified Lambert-Beer equation: 



N x = 



lnlO jAdv 



c ret 
^X 



(1) 



where the lnlO is needed to convert the integrated absorbance, A, 
to optical depth and S^' is the experimental RAIR band strength 
of species X. Transmission band strengths available from the lit- 
erature cannot be used, because the total number of molecules 
probed with RAIRS cannot directly be related to a value probed 
in a transmission absorption experiment. This is because the in- 
cident beam goes through the ice layer twice with an angle to 
the surface. Instead, values of S J f have been calculated from a 
calibration experiment without hydrogenation where the depo- 
sition rate is 1.0xl0~ 7 mbar s _1 to grow a layer of typically 10 
to 20 ML and where the sticking probability is assumed to be 
1. It should be noted that these experimental values can differ 
between different experimental set-ups and even for the same 
set-up over time, because they are determined by intrinsic prop- 
erties such as the alignment of the system. They should therefore 
not be used to compare directly with observations of interstel- 
lar ices and even a comparison between different experiments 
should be made with caution. The values for S^ f as well as the 
spectral assignments of the vibrational modes are summarized 
for all species in Table [2] 
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Table 2. Overview of the integrated frequency ranges with cor- 
responding spectral assignment and band strength. The error on 
S r v ef amounts to -30-40%. 



Species 


Mode 


Integration range 


c ret 
°X 






(cm- 1 ) 


(cm molecule -1 ) 




v 3 


2265-2245 


1.3(-17) 


13 c 18 o 2 * 




2265-2245 


8.8C-17) 


13 C 18 0/ 




2265-2245 


4.8(-17) 


13 C 18 02 rf 




2265-2245 


4.3(-17) 


13 C 18 0/ 




2265-2245 


3.2(-17) 


13 C 18 02 / 




2265-2245 


2.7(-17) 


CO" 


Vl 


2160-2120 


9.5(-18) 


co' 




2160-2120 


6.4(-18) 


HCOOH 


V3 


1800-1550 


1.3(-16) 


CH 3 CHO 


v 7 


1745-1719 


8.0(-18) 


CH^CHO 




1365-1330 


2.8(-18) 


CH,OH 


vs 


1060-990 


l.K-17) 


CH 4 




1320-1290 


8.2(-18) 



"pure 13 C 18 2 , "22:78% 13 C 18 2 :H 2 0, '39:61% lj C 18 2 :H 2 0, 
rf 48:52% 13 C 18 2 :H 2 0, "45:55% 13 C 18 2 :CO, -^80:20% 13 C 18 2 :CO. 



There are several contributions to the uncertainty in the 
experimentally measured band strengths. The largest fraction 
comes from the actual deposition and sticking onto the surface. 
Since our experiment does not contain a micro-balance, the ab- 
solute number of molecules on the surface is not known and the 
possibility that some molecules freeze out on other surfaces than 
the Au substrate cannot be ruled out. This leads to systematic 
errors for but there is no effect on the relative error for 
between different experiments with the same ice morphology. 
Based on experiments aimed at the same ice thickness, inaccu- 
racies in the deposition flow are estimated to be ~30%. Since 
the integrated area of the absorption band is determined very 
precisely for the deposited species, the column densities can 
be accurately normalized, but the derived band strengths can- 
not. Another source of uncertainty for ice mixtures is the preci- 
sion of the mixing ratio, which is of the o rder of 10% (see also 
Qber g et ail 120071: [Bisschop et al.Ll2007al) . In summary, the ac- 
tual uncertainty on the band strengths is substantial and ranges 
from 30^40%, but the relative uncertainty for ices with the same 
morphology is less than 5%. 

3.2. Reaction rate calculations 

Th e method fo r calc ulating reaction rates is described in detail 
by iFuchs et al.l d2007l) . In short, a species X can react with H- 
atoms to form species Z through: X + H i", Z. The column den- 
sity of X that has reacted, 2Vx(0> is given by: 



dN x (f) 
dt 



-k NnN x , 



(2) 



where ko stands for the reaction constant and Nn for the sur- 
face density of H-atoms. In our experiment the H-atom flux is 
kept constant and corresponds to 7.8xl0 13 or 5.0xl0 14 cm -2 s _1 . 
Furthermore, the atoms have a certain penetration depth 
as observed in experiments of H-atom reactions with CO 
dWatanabe et al.L 12004 IFuchs et aTll2007l) which means that not 
all of the deposited parent species are available for the reaction. 
Consequently, Nx(t) is calculated via: 



where ao is the fraction of Nx(0) that is available to react and 
fio (in min -1 ) corresponds to koNn/60 (the factor 60 comes from 
the conversion of seconds to minutes). In specific cases only up- 
per limits on reaction rates are calculated and then it is assumed 
that cq only includes the uppermost ice layer. Limits based on 
the column density decrease after 1 minute give the most con- 
servative upper lim it on fio- Similar to the cases described by 
IFuchs et aU d2007l) . fits to the reaction rate differ when they are 
made over the complete time-period of the experiment (hrs) or 
only over a shorter period (minutes). 

3.3. TPD analysis and calculation of the production yield 

The TPD data provide complementary information on the reac- 
tion yields and are important, in particular, for those molecules 
that are not accurately determined by RAIRS. The TPD data 
are fitted by second order polynomial baselines. The tempera- 
ture range over which a baseline is fitted depends on the desorp- 
tion temperature of a specific species. Calibration experiments 
of pure ices with a known number of molecules have been per- 
formed by measuring the corresponding integrated area of the 
mass spectrometer signal. The number of molecules Nz in other 
experiments has been determined by comparison of the inte- 
grated signal to the number of molecules in the calibration ex- 
periments. Since pure ices are needed for the calibration, no ac- 
curate yields can be calculated for H2CO which is not readily 
available as a pure ice due to polymerization (see § 16.11 ). The 
yield, Yz, of the newly formed species Z in % can then be calcu- 
lated through: 



N z 



aoiVx(O) 



(4) 



N x (t) = N X (Q) a (1 - e**) 



(3) 



where Nz is the column density as derived from the TPD data, 
ao is taken from the fit of the RAIRS data and Nx(0) is the ini- 
tial column density of the precursor species as derived from the 
RAIR spectra. 

4. CO? containing ices 

4.1. Results 

Pure ices of CO2 are bombarded with H-atoms to search for the 
possible formation of HCOOH following the reaction route as 
shown in Fig.Q] The reactivity of CO2 can be tested by record- 
ing a decrease in the CO2 RAIR signal and by monitoring any 
reaction products. Unfortunately, the CO2 2300 cm -1 ice band is 
difficult to quantify, because this band overlaps with rotation- 
vibration transitions of CO2 present in the purge gas in our 
spectrometer. The use of I3 C 18 02 isotopic species does not im- 
prove this situation. Therefore we have focused on the strongest 
HCOOH band, the C=0 stretching mode, at 1710 cm -1 , to mon- 
itor formic acid formation. The observed difference RAIR spec- 
trum of 13 C O2 bombarded by H-atoms has been compared to 
that of 12 C 18 02 in Fig. [2^. Both spectra show very weak fea- 
tures at 1730 cm -1 and 1500 cm -1 , but none at positions typical 
for HCOOH. The detected bands, however, occur at exactly the 
same positions as previou sly seen for H2CO whe n formed upon 
CO hydrogenation (e.g., Wat anabe et al.L 120041) . Furthermore, 
the features do not shift when a different isotopic species is used, 
which is expected when the formation involves CO2 ice. In the 
TPD spectra of CO2 ices bombarded with H-atoms for 3 hrs (not 
shown) mass 29 and 30 amu desorb in two steps around 100 K 
and 140 K, and the formation of other species is not detected 
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Table 3. Upper limits on the reaction/destruction rates for 
HCOOH and CO2. The uncertainty on Icq amounts to a factor 
2. 



Species 


Ice matrix 


T 








(K) 


(cm 2 s~') 


co 2 


pure 


14.5 


<6.2(-17) 


C0 2 


13 C 18 2 :H 2 39:61% 


14.5 


<6.0(-17) 


C0 2 " 


13 C 18 2 :H 2 39:61% 


14.5 


<3.8(-17) 


C0 2 


13 C 18 2 :H 2 22:78% 


14.5 


<6.7(-17) 


C0 2 


13 C 18 2 :H 2 48:52% 


14.5 


<3.2(-17) 


C0 2 
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" Limit derived for experiment without bombardment. 



through TPD. These t emperatures and m asses are identical to 
what was observed by Fuch s et ail d2007l) for H2CO desorption 
from CO ices bombarded by H-atoms. During H-atom bombard- 
ment an increase in the mass 28 amu signal is observed, which 
is likely due to degassing of both CO and N2 from the metal 
parts of our experiment. It is therefore plausible that the mea- 
sured low level of H2CO formation observed in the experiment 
originates from hydrogenation of background gaseous CO and 
is not related to the CO2 ice. 

To test whether the presence of H2O affects the reactivity 
of CO2 ice upon H-atom bo mbardment as has been observed 
for CO in CO:H 2 mixtures dWatanabe et aUl2004lFuchs et aU 
120071) , mixtures of 22:78% to 48:52% C0 2 :H 2 have been in- 
vestigated (see Fig. Wp). Like in the experiments with pure CO2, 
weak RAIR features of similar intensity are observed at 1730 
and 1500 cm . Again we assign these features to the C=0 
stretching mode and the C-H bending mode of H2CO. Thus, as 
for the pure ices, a small amount of background CO accretes and 
forms H2CO. Within the sensitivity of our experiment, we con- 
clude that CO2 does not react with H-atoms even if mixed with 
H 2 0. 

Finally, mixtures of 12 C 16 and 13 C 18 02 are studied to de- 
termine which species is more likely to react upon H-atom 
bombardment: CO or CO2. Since CO h ydrogenation reactions 
were previously re ported in the literature (Wata nabe et alll2004t 
Fuch s et all 120071) . the answer to this question must be CO. In 
Fig. [2J: and d, the resulting difference spectra for C02:CO mix- 
tures are shown for different ice thicknesses and with mixture 
concentrations of 45:55% C0 2 :CO and 80:20% C0 2 :CO, re- 
spectively. Similar to H-atom bombardment of pure CO2 ices 
and C02:H20 mixtures, no evidence for HCOOH formation is 
observed in C02:CO mixtures. In contrast, CO does react with 
H-atoms to form H2CO and CH3OH as is evidenced by the pres- 
ence of strong H2CO absorption features at 1730 and 1500 cirT 1 
and CH3OH at 1030 cm -1 . This is consistent with H-atom bom- 
bardment experiment s for pure CO and CO :H2Q mixtures by 
Watanabe et alj d2004l) and lFuchs et alj d2007l) . The complemen- 
tary TPD data show the same picture of no HCOOH formation 
and clear H2CO and CH3OH formation from CO. Other prod- 
ucts than the precursor and product species H2CO and CH3OH 
are not observed. 

4.2. Reaction rates 

No hydrogenation products of CO2, specifically HCOOH, 
are observed within the experimental sensitivity. The limit 
on the formation reaction rate for HCOOH from CO2 is 
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Fig. 2. Difference spectra of the 1800-1400 cirT 1 range after 
180 minutes of H-atom bombardment, (a) 15 ML pure C 18 02 
ice (black line) and 13 C 18 2 (grey line), (b) 15 ML C0 2 :H 2 
39:61% (black line), 22:78% (grey line) and 48:52% (light grey), 
(c) 45:55% C0 2 :CO mixtures for 15, 30 and 45 ML ice thick- 
ness and (d) 80:20% C0 2 :CO mixture for 15 and 30 ML. The 
temperature of the ice is ~14.5 K. The arrows indicate how the 
absorbance of the H2CO bands decreases with decreasing thick- 
ness. 
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Fig. 3. A/A for the CO 2140 cm" 1 band in C0 2 :CO mixtures at 
14.5 K for the 45:55% (upper panel) and 80:20% (lower panel). 
The symbols refer to 15 ML +, 30 ML *, and 45 L O. The dotted 
lines indicate the fits to the data. 



<7.0xl0~ 17 cm 2 s based on the limit on the column density 
for HCOOH after 1 min of H-atom bombardment for all ice 
morphologies (see § 13.21 for the derivation and Table [3] for the 
individual values for each experiment). 

In Figure|3]the absorbance divided by the initial absorbance 
at t =0, A/Aq, is shown for the 45:55% and 80:20% C0 2 :CO ice 
mixtures. The data are fitted as described in § |3.2| and the fits are 
indicated in Fig.[3]with dotted lines. The resulting values for ao 
and /3q as well as the ko are given in Table [4] Since the H2CO 
band strength could not be determined accurately in these exper- 
iments only ao and ySo are fitted. Clearly ao, i-e., the fraction of 
CO molecules available for reaction, decreases with increasing 
ice thickness. 



6 



Bisschop et al.: H-atom bombardment of COt, HCOOH and CH3CHO ices 



Table 4. Values for oro, fio, and the reaction rate £0 for CO in 
CO:C02 mixtures at ice temperatures of 14.5 K upon H-atom 
bombardment. The uncertainties for o-q and /3q amount to 10- 
20% and for ko are a factor 2. 



Ice mixture 


Ice thickness 


a 


A) 


k 




C0 2 /total 










(ML) 




(min" 1 ) 


(cm 2 s" 1 ) 


45:55% C0 2 :CO 


8/15 


0.75 


0.024 


2.9(-15) 


45:55% C0 2 :CO 


17/30 


0.51 


0.018 


2.2(-15) 


45:55% C0 2 :CO 


25/45 


0.34 


0.014 


1.7(-15) 


80:20% CO: CO 


3/15 


0.67 


0.032 


3.8(-15) 


80:20% CO.: CO 


6/30 


0.57 


0.017 


2.0(-15) 



4.3. Discussion and conclusion 

Previously, CO2 + H r eactions have been studied by 
Milligan & Jacox| (Il97ll) in UV-photolysis experiments of 
Ar:C02:H20 matrices. Only the formation of more oxygen-rich 
species such as CO3 has been observed, but not HCOOH. 
Although these results cannot be compared directly with ours 
(see § |2.2| >, our findings are consistent with theirs that CO2 
does not react readily with H-atoms. To explain the lack of CO2 
hydrogenation reactions, Fig. |4] presents the relative formation 
energies of possible products. Reaction of CO2 with H-atoms to 
either HCO+O or CO+OH is energetically highly unfavorable. 
This is due to the HOCO transition state being ~ 15930 K 
1 .4 eV) higher in energy compared to CO2+H in the gas phase 
Lakin et all 120031) . On the other hand, a hydrogenation reaction 
could be expected based on the higher heat of formation of 
C0 2 +2H with respect to HCOOH of -49370 K (4.3 eV). 
Indeed, in the chemical physics literature CO2 is found to 
hydrogenate to HCOOH on ruthenium and iridium catalyst s 
through formate complexes with the catalyst dOgo et all 12006). 
The chemically bonded formate species subsequently reacts 
with H30 + to form HCOOH. This reaction mechanism requires 
acidic species as well as a catalytic surface that ar e not present 
in the current experiment. Hwang & Mebel d2004l) calculated a 
potential energy surface for the gas phase H2 + CO2 reaction. 
The HCOOH end product is higher in energy than the initial 
species by -2600 K (0.22 eV), but the C0 2 + 2H reaction is 
exothermic (see Fig . |4j) . For the reaction of CO2 with H-atoms or 
H2 the same transition state H2CO2, a complex cyclic structure, 
has to be overcome. This transition state lies 35000-37000 K 
(3.0-3.2 eV) above the starting point, CO2 + H2, but below 
CO2 + 2 H. This step can therefore not be rate-limiting for the 
CO2+ 2H — » HCOOH reaction. Since the reaction is clearly not 
observed in our experiment, the rate determining step must be 
the addition of the first hydrogen atom to CO2 to form HCO2 
and the barrier to this reaction must be too high to be overcome 
for the ice temperatures of 12-60 K as used in our experiments. 

Several recent studies have focused on reactio ns of CO + 
H lead ing to the formation of H2CO and CH3OH. iFuchs et all 
(2007) find that even for the lowest ice thicknesses of 1 to 2 ML 
30% of the ice is hidden from the impinging H-atoms. At higher 
thicknesses ao increases and there is a maximum layer thickness 
of 12 ML of CO ice that can react with H-atoms. The behavior 
of cko for our C02:CO mixtures is consistent with this picture, 
although at maximum 7+3 ML and 8+3 ML of the CO ice reacts 
for the 45:55% and 80:20% C0 2 :CO mixtures, respectively. In 
other words mixing CO with CO2 does not cause more CO to be 
"hidden" from the H-atom exposure. 

Our reaction rates, Pq, of e.g., 0.032 min" 1 in CO:C02 
80:20% with 15 ML total ice thickness at 14.5 K are simi- 
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Fig. 5. The difference spectrum of the HCOOH v s (C=0) stretch 
for HCOOH bombarded for 4 hrs with H-atoms at 12 K (solid) 
and at 40 K (dashed). 

lar to those found for pure CO ices by IFuchs et al.1 d2007l) of 
~0.030 min" 1 (~0.023 mkr 1 when converted to our assumption 
for the initial number molecules) for a similar amount of CO of 
1 1 ML at 15 K. These rates are the same wit hin the 30-40% un- 
certainty. However, forCO:H20 1:5 mixtures lFuchs et ail d2007l) 
find a value of -0.1 1 min" 1 (0.083 min" 1 ) for 12 ML at 15 K 
indicating that the reaction rates for CO hydrogenation are sig- 
nificantly higher in mixtures with H2O ice. The similarity be- 
tween the reaction rate of CO in mixtures with CO2 and pure 
CO ices and the difference between those and CO:H20 ice mix- 
tures can be explained by CO and CO2 only interacting through 
weak Van der Waals forces. The electronic structure of the CO 
molecule will therefore not differ significantly in mixtures with 
CO2 from pure CO ices. H2O on the other hand has a stronger 
dipole moment of 1 .85 D compared to zero and 0. 1 1 D for CO2 
and CO, respectively, and forms hydrogen bonds. Furthermore 
it is known that CO strongly int eracts with and influen ces the 
band strengths of H2O molecule (Bouwman et all [2007b . Thus, 
the electronic structure of the CO molecule will be perturbed in 
mixtures with H2O, strongly affecting the reaction rate of CO 
with H-atoms. In summary the presence of CO2 in ice mixtures 
with CO does not strongly affect the reactivity of CO with H- 
atoms. 



5. HCOOH containing ices 

5.1. Results 

Figure [5] shows the difference spectrum for vs(C=0) at 
~1710 cm" 1 of pure HCOOH ice bombarded with H-atoms as 
well as control experiments with bombardment of H2 molecules 
at 1 2 K (for an overview of all infrared features of HCOOH 
see [Cyriac & Pradeepl 120051) . The growth of an infrared fea- 
ture around ~1050 cm 1 indicative for CH2(OH)2 formation 
has not been obser ved (for an overvie w of the infrared features 
of CH 2 (OH) 2 see iLugez et all 1 19941) . A decrease on the blue 
side of the v s (C=0) mode at 1710 cnr 1 of HCOOH is seen at 
1750 cm" 1 as well as an increase at 1730 cm" 1 , which means that 
the overall HCOOH band profile changes slightly. At 1730 cm" 1 , 
the C=0 stretch for H2CO is located, but other features of H2CO, 
such as the 1500 cm" 1 band, are missing. The decrease cor- 
responds to <0. 1 ML derived from our calculated RAIR band 
strength. These features are present in difference spectra for 
HCOOH ice bombarded with H-atoms at 12 and 40 K. A simi- 
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Fig. 4. Potential energy scheme for CO2 dissociation and hydrogenation. The relative energies are b ased on the heat s of formation 
at K. The scal e in electronyolts is indicated on the right. The heats o f formation are derived from ICox etalJ(ll989h for H . C0 2 , 
CO and O, from lRuscic eTaT] d2002l) for OH, from lGurvich et all dl989b for HCO and HCOOH and from lLakin et alj d2003b for the 
HOCO transition state(ts). 



lar shift is seen for transmission infrare d experiments with pur e 
HCOOH ice that is heated to -60 K dBisschop et all l2007ah . 
At the same time the vg(CH) and vs(OH) vibrational modes in- 
crease due to conversion of HCOOH in dimeric form to HCOOH 
organized in chains. In the RAIRS spectra these bands are also 
seen to increase. Furthermore, the same change in RAIR profile 
is found for HCOOH ices of 40 K, where H-atoms cannot stick 
any longer onto the surface, but can only collide. Since the ice 
has a temperature of 40 K, the reorganization of the ice is less 
and consequently the signal of the difference spectrum is smaller. 
In conclusion, the RAIR data do suggest that some restructuring 
takes place in the surface but no reaction. 

With TPD the masses of 48 (CH 2 (OH) 2 ), 46 (HCOOH), 45 
(HCOO), 44 (C0 2 ), 32/3 1 (CH 3 OH), 30/29 (H 2 CO), and 28 amu 
(CO) have been monitored during warm-up. No products are de- 
tected at 48, 32, 31, or 30 amu to upper limits of <0.01 ML, 
indicating that HCOOH is neither hydrogenated nor dissociated. 
Thus, consistent with the lack of a 1500 cirT 1 H 2 CO absorption 
feature in the RAIRS data, no evidence for H 2 CO formation is 
observed in the TPD experiment. The detected masses 45, 44, 
and 29 amu are assigned to HCOOH dissociating in the mass 
spectrometer, because the same relative mass ratios are seen for 
a TPD spectrum of pure HCOOH ice that is not bombarded by 
H-atoms. We conclude that within the limits of our experimental 
set-up the reaction of HCOOH with H-atoms is not efficient at 
12 K. 

5.2. Reaction rates 

Since no unambiguous evidence for HCOOH destruction in the 
ice is found, it is only possible to derive an upper limit on its 
reaction rate, presented in Table [3] It is clear that the HCOOH 
destruction rates are below 2.3xl0~ 17 cm 2 s _1 as derived from 
the limit on the column density after 1 min of H-atom bombard- 
ment (see § 13.21 ). As for C0 2 these reaction rates are very low. 



5.3. Discussion and conclusion 

In chemical physics literature HCOOH hydrogenation on cat- 
alytic surfaces has been shown to lead to decomposition 
of HC OOH rather than methanediol formation dBenitez et all 
119931) . HCOOH adsorbs onto such a surface as HCOO" and H + 
which can be further hydrogenated. The catalytic surface, how- 
ever, clearly affects the end products and overcomes a reaction 
barrier that prohibits spontaneous decomposition. If hydrogen 
atom addition and dissociation occur simultaneously, C-0 bond 



10480 K H-CO + OH + H 
HCOOH + 2H HCO + aO + H 



5090 K 



0.0 eV 
-0.4 eV 



HfO + H,0 



Fig. 6. Potential energy scheme for HCOOH hydrogenation. The 
energies are based on the heats of formation at K. The en- 
ergy scale in electronvolts is i ndicated on th e right . The heats 
of formation are derived from Gurvich et al. (1989) for HCO, 
H 2 CO and HCO OH, from lRuscic et all (120021) for OH and from 
ICoxetalJ (119891) for H and H 2 0. 



cleavage is more energetically favorable (as shown in Fig. [6j. 
However, it is clear from the results in § 15.11 that no H 2 and 
H 2 CO formation occurs. Thus a high barrier for H-addition to 
HCOOH must exist for both mechanisms and HCOOH + H re- 
actions in the ice are inefficient. 



6. CH3CHO containing ices 

6.1. Results 

The infrared spectroscopic features detected for pu r e CH3 CHO 
ice match with those detecte d by iBennett et all (I2005I) and 
iMoore & Hudsor3dl998il200l . The strongest CH 3 CHO band is 
the C=0 stretching mode, v s (C=0), at 1728 cirT 1 (5.79 yum). 
During H-atom bombardment the intensity decreases, but a 
small positive wing is observed at 1710 cirT 1 (see Fig.|7]i. This 
band is assigned to the C=0 stretching mode of H 2 CO. Other 
CH3CHO features (e.g., the umbrella deformation mode, vd, at 
1345 cm 1 ) also decrease and new bands appear at 1030 and 
1300 cirT 1 that are attributed to the C-O stretching mode of 
CH3OH and the deformation mode of CH4, respectively. No 
clear absorption is observed at 1050 cm -1 , where the strongest 
C 2 H50H band, the C-O stretching mode, is expected. Since 
this frequency region is particularly problematic in our detec- 
tor, the detection upper limit on A^(C 2 H50H) amounts to only 
3xl0 15 molecules crrT 2 , i.e., 3 ML. Another strong band of 
C 2 H50H is expected at 3.5 fj.m. Unfortunately, this feature over- 
laps with a number of CH3OH modes. Broad weak features are 
indeed detected in this range, but due to the complexity of both 
C 2 H50H and CH3OH absorptions and the relatively weak signal 
this cannot be used to determine whether C 2 H50H is present. 
Additionally, it is important to note that no strong features are 
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Fig. 7. The difference spectrum of CH3CHO after 3 hrs of 
bombardment with H-atoms. The negative peaks correspond to 
CH3CHO destruction, the positive wing around 1700 cm 1 is 
assigned to the vs(C=0) mode of H2CO and the features at 
1300 cm 1 and 1030 cm 1 to v D (CH 4 ) of CH 4 and v s (C-0) of 
CH3OH, respectively. The arrow indicates the position where the 
strongest C2H5OH absorption is expected. The dotted line in- 
dicates the RAIR spectrum for pure CH3OH, showing that the 
feature detected at 1030 cm 1 matches that of pure CH3OH. 
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Fig. 8. The TPD spectrum for 40 ML of CH3CHO bombarded 
with H-atoms for 3 hrs at 14.5 K. The black line refers to the 
16 amu signal (CH4), dark grey line to 31 amu (CH3OH), grey 
to 46 amu (C 2 H 5 OH) and light grey to 44 amu (CH3CHO). 



observed around 2140 cm , where both CO and CH2CO have 
infrared features. This is perhaps not surprising, because the for- 
mation of CO would involve not only the breaking of a C-C 
bond, but also hydrogen-abstraction, which is not very likely 
in this hydrogen-rich environment. The formation of ketene, 
CH2CO, is even less likely because its formation is strongly en- 
dothermic. 

The formation of CH4, H2CO and CH3OH is corroborated 
by the TPD spectra, where 16, 30 and 31 amu mass peaks at 
45 K, 100 K and 140 K are found, respectively (see Fig. [8] for 
CH4 and CH3OH). The peaks for 16 amu at higher temperatures 
are due to O-atoms detected by the mass spectrometer when 
other molecules dissociate. The desorption temperatures for 16 
and 31 amu are similar to those of pure CH4 and CH3OH ice 
confirms their RAIR detection. The TPD spectra and desorption 
temperatures of 29 amu are co nsistent with the desorp tion tem- 
peratures for H 2 CO found by IWatanabe et all d2004l) . In addi- 



Fig.9. The A/A ratio (upper panel) and AA(CH 3 CHO) (lower 
panel) for the CH3CHO 1345 cirT 1 band for different ice thick- 
nesses and a constant ice temperature of 14.5 K. The symbols 
refer to 11.4 ML (+), 7.8 L (*), 21.2 L (O), 45.8 ML (□) and 
56.0 ML (a). The dotted lines indicate the fits to the data. 
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Fig. 10. The A/A ratio (upper panel) and AJV(CH 3 CHO) (lower 
panel) for the CH3CHO 1345 crrT 1 band for different ice tem- 
peratures and a constant ice thickness of 1 1.3 ML. The symbols 
refer to 14.5 K (+), 12.4 K (*), 15.8 K (□), 17.4 K (O) and 19.3 K 
(a). The dotted lines indicate the fits to the data. 



tion, a TPD desorption peak is located at ~160 K for masses 45 
and 46 amu (see Fig . [8j ■ This is assigned to C2H5OH desorption 
based on a comparison with the TPD of pure non-bombarded 
C2H5OH ices. In summary, a fraction of CH3CHO, below the 
infrared detection limit of the 1050 cirT 1 band, is converted to 
C2H5OH and a larger fraction forms CH 4 , H 2 CO and CH3OH. 
So even though the conversion of acetaldehyde to ethanol is not 
complete, it is important t o note that a pathway in the proposed 
hydrogenation scheme by iTielens & Charnl evld 19971) is experi- 
mentally confirmed. 

6.2. Reaction rates and production yields 

The value for A^CFLCHO) as derived from the v D (umbrella) 
spectral feature at 1345 cirT 1 is shown in Fig. [9] as a function 
of time for different ice thicknesses at 14.5 K. Also shown are 
the fits to the data. The VD(umbrella) mode is chosen for analy- 
sis rather than the 1728 cirT 1 band, because the latter overlaps 
with the v s (C=0) of H 2 CO at 1720 cm -1 . Clearly, the absolute 
amount of CH3CHO that can react increases with ice thickness, 
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Fig. 11. The ao and /3o dependencies of the CH3CHO+H reaction 
on thickness and temperature. The temperature of the ices for the 
thickness dependence experiments is constant at 14.5 K. The ice 
thickness is similar for the temperature dependence experiments 
at 11-12 ML. 

whereas A/Aq decreases. The temperature behavior is more com- 
plex and is shown in Fig.fTOl 

The ffo and fto values derived from the fits as function of the 
thickness and temperature are shown in Fig. [TT] The values for 
ao decrease with increasing thickness, but do not depend on ice 
temperature within the measured regime. The latter is not sur- 
prising as the CH3CHO ice structure does not change between 
1 5 and 75 K. The value for flo is independent of ice thickness, but 
does depend on ice temperature. It is largest for ice temperature s 
between 15-16 K, similar to the case of CO dFuchs et all 12007). 
This is expected as the maximum reactivity is mostly determined 
by the mobility of H-atoms at the surface. At low temperatures 
H-atoms move more slowly resulting in a lower reaction rate. At 
higher temperatures the diffusion rate is higher but has to com- 
pete with an increased evaporation rate. 

For the C2H5OH formation, only yields can be calculated 
from the TPD data because the RAIR feature at 1050 cm 1 
overlaps with the v s (CO) band of CH3OH. The yields for CH 4 , 
CH3OH and C2H5OH are given in Table [5] Even when consid- 
ering that there is a general quantitative uncertainty of ~10%, 
it is clear that the summed yield of the different products is 
not 100%. This is most likely due to missing H2CO yields, be- 
cause these are not reliably calibrated. Furthermore IXCFU) is 
expected to be equal to F(H2CO+CH 3 OH), because CH3OH is 
formed from H2CO after CH3CHO dissociation. However, the 
CH3OH yield is significantly higher than CH4 (see § I6.31 >. The 
solid state C 2 H 5 OH yields are <20%. 

6.3. Discussion and conclusion 

Previously, C2H5OH and CH3CHO were shown to form in in- 
terstellar ice analogues by photoly sis of C2H2:H20 mixtures 
dMoore et all 1200 UlWu et alll2002l) . Since in such experiments 
both OH and H fragments are present with excess energy, it is 
difficult to disentangle potential pure hydro genation reactions 
and re actions involving OH radicals. Indeed. lMoore & Hudson! 
(2005) explain formation of C 2 H 5 OH and CH3CHO by reactions 
of C2H5 and C2H3 with OH, respectively. In this paper we focus 
on the reactions with thermal H-atoms only. 

Since H2CO is known to react with H-atoms to CH3OH 
dWatanabe et all |2004 iHidaka etall |2004|) it is thus likely that 
the next more complex aldehyde, acetaldehyde (CH3CHO), will 



form ethanol (C2H5OH). In Figure Q~2] the relative heats of 
formation at K are show n dWiberg et all [1991 ICox et all 
1989; iGurvich et all 119891: iFrenkel et all Il994t iMatus et al. . 
2007). The exothermicity of CH3CHOH formation is higher 
than that of CH3CH2O, but which of the species is more likely 
formed depends on the reaction barriers. Subsequent formation 
of C2H5OH is likely fast, because reactions of radicals with H- 
atoms commonly have no activation barriers. As described in 
§ l6.2l onlv a fraction of CH3CHO is converted to C2H5OH, and 
a larger fraction leads to CH4, H2CO, and CH3OH formation. 
For hydrogenation of CH3CHO a C=0 bond is converted to a 
C-O bond instead of breaking a C-C bond. Since the C=0 bond 
is intrinsically stronger it is likely that the entrance channel to 
hydrogenation is higher in energy compared to dissociation. 

Thus H-atoms can break the C-C bond as well as the C=0 
bond to form CH 4 , H 2 CO and CH 3 OH or C 2 H 5 OH in ices as 
prepared here. As shown in Fig.[T2]the formation of CH4+HCO 
is more exothermic than that for CH3+H2CO. Furthermore, the 
energy released in this step is higher than the binding energy 
of CH 4 to the surface, which is -700 K (0.06 eV). This likely 
explains why the T(CH 4 ) is lower than F(CH 3 OH+H 2 CO), be- 
cause the formation energy is sufficient for CH4 desorption. The 
energy released during the formation of H2CO, CH3OH and 
C2H5OH is even higher and may also cause a fraction of the 
molecules to desorb. 

7. Astrophysical implications 

Our experiments show that CO2 reaction rates with H-atoms are 
very low, making it an implausible route for HCOOH forma- 
tion. A numbe r of other HCOOH formation routes are possi- 
blejseee^. iMilligan & Jacoxl I197U iHudson & Moorel 119991: 
lKeanell200ll) . from either HCO+OH -> HCOOH or HCO+O -> 
HCOO+H — > HCOOH. In addition, experiments suggest that un- 
der specific catalytic conditions CO2 can react to form HCOOH 
dOgo et all 120061) but this requires catalytic surface sites, i.e., 
CO2 directly attached to a silicate or metallic grain site. Such 
a situation is less likely in dense interstellar clouds where thick 
ice layers have already formed and cover any potential catalytic 
sites. In conclusion, under astrophysically relevant conditions 
solid CO2 in bulk ice is a very stable molecule that is not likely 
to react with H-atoms. 

Similar to CO2, reaction rates of HCOOH with H-atoms are 
below the detection limit in our experiment. Formation of the 
so far undetected interstellar species CH2(OH)2 in this way thus 
seems unlikely. Unless other formation mechanisms are found 
an observational search for this species based upon solid state as- 
trochemical arguments is not warranted. We conclude that CO2, 
HCOOH and CH 2 (OH) 2 do not appear to be related through 
successive hydrogenation in interstellar ice analogues under the 
conditions as used in the present study. 

In contrast to C0 2 and HCOOH, CO does react with H- 
atoms. The reaction rates of CO in CO:CQ2 mixtures are very 
similar to those found by iFuchs et ail d2007l) for pure CO ices. 
CO hydrogenation in interstellar ices will thus not be strongly 
affected by the presence of CO2 in the ice. It is likely that the 
reaction rate is the same for H+CO independent of the CO con- 
centration and that of other species in apolar interstellar ices. 

Reactions of CH3CHO in interstellar ices proceed at simi- 
lar rates compared to CO hydrogenation. A maximum of 20% 
will be converted to ethanol, C2H5OH, while another major re- 
action channel leads to CH4, H2CO and CH3OH The precise 
abundance of CH3CHO in interstellar ice is not yet well deter- 
mined. However, abundances of 1-5% are quoted in the litera- 
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Table 5. Values for chq, (3q, the reaction rate Icq for CH3CHO, and production yields, Y(X) upon H-atom bombardment. The uncer- 
tainties for ao and /3o amount to 10-20%, for Icq are a factor 2 and Y 20%. 



Ice temperature 


Ice thickness 




A> 


ko 


K(C 2 H 5 OH) 


y(CH 4 ) 


y(CH 3 OH) 


(K) 


(ML) 




(min ') 


(cm 2 s~') 


(%) 


(%) 


(%) 


14.5 


7.8 


0.40 


2.3(-2) 


2.8(-15) 


15 


21 




14.5 


11.4 


0.36 


1.5(-3) 


1.8(-16) 


10 


19 




14.6 


13.5 


0.31 


2.2(-2) 


2.6(-15) 


14 


22 


48 


14.5 


18.8 


0.30 


1.9(-2) 


2.3(-15) 


13 


18 




14.5 


21.2 


0.18 


2.1 (-2) 


2.5(-15) 


17 


24 




14.6 


22.1 


0.22 


2.0(-2) 


2.4(-15) 


20 


23 


39 


14.6 


45.8 


0.13 


1.8(-2) 


2.2(-15) 


20 


21 


38 


14.6 


56.0 


0.13 


1.7(-2) 


2.0(-15) 


20 


17 


35 


15.8 


11.6 


0.39 


1.6(-3) 


1.9(-16) 


21 


22 


33 


17.4 


11.3 


0.45 


9.9(-3) 


1.2(-15) 


16 


19 




19.3 


11.2 


0.39 


1.1 (-2) 


1.3(-15) 


14 


13 


15 



Fragmentation: 
CH 3 CHO + 2H U /U \ Hj+ ^ co + H 



9720 K 



CH 4 +HCO + H 

43600 K 



Hydrogenation: 

6310 K 

CH 3 CH 2 + H 



CH 4 + H,CO 



CH3CHO + 2H 



11650K 



CFLCHOH + H 



47000 K 



CHjOH 



. 0.0 eV 

- -0.5 eV 

- -1.0 eV 



-5.1 eV 



Fig. 12. Potential energy scheme for CH3CHO fractionation and hydrogenation. The relative energies are based on the heats of 
formations at K. An approximat e energy scale in electro nvolts is given on the right. The heats of formation are derived from 
IWibere et all (fl99ll) for CH3CHO, ICoxet al.1 (fl989h for H, iGurvich et all (fT989b for CH 3 , H 2 CO, CH4 and HCO, iFrenkel et alJ 
(11994 for C 2 H 5 OH and lMatus etaD (l2007h for the CH3CHOH and CH3CH2O radicals. 



ture ( 


Schutte et all 1 1 997L 1 1 999b iGibb et al.Ll2004tlBooeert et all 


|2004 


). These values can be used to derive an upper limit on 



the C2H5OH abundance that could thus be formed. As an ex- 
ample we compare the abundances for the high mass source 
W 33A, where CH3CHO has a solid state abundance of 9.8xl0~ 6 
and CH3OH of 1.4-1.7xl0~ 5 both with respect to H 2 . If we 
assume that all solid CH3CHO is present in the surface layer 
and the C 2 H 5 OH yield is -20%, an abundance of C 2 H 5 OH 
of at most 2.0xl0~ 6 with respect to H 2 can be formed. This 
leads to an upper limit on the C 2 H 5 OH/CH 3 OH ratio of 0.14. 
In reality this value will be lower as part of the CH3CHO ice 
may be shielded from incoming H-atoms and other destruc- 
tion reactions will likely be competing with hydrogenation re- 
actions. The limit of 0.14 is clearly higher than the observation- 
ally derived C 2 HsOH/CH30H abundan ce ratio in the gas phase 
of 0.025+0.013 dBisschop et al.L l2007b). Formation of C 2 H 5 OH 
from solid state hydrogenation of CH3CHO is thus sufficient to 
explain the observed abundances of C 2 H50H. 

8. Summary and conclusions 

Hydrogenation reactions of C0 2 , HCOOH and CH3CHO inter- 
stellar ice analogues have been studied under ultra-high vacuum 
conditions. RAIRS and TPD have been used to analyze the re- 
sults. From these experiments reaction rates and upper limits on 
destruction and formation rates of the above mentioned species 
are calculated. The main conclusions derived from this work are: 

- C0 2 and HCOOH do not react with H-atoms at a detectable 
level. Only minor fractions of the species desorb due to the 
bombardment. Solid state formation of HCOOH from C0 2 
and CH 2 (OH) 2 from HCOOH are likely inefficient in inter- 
stellar ices. 



- Hydrogenation of CO to H 2 CO and CH3OH from CO mixed 
with C0 2 has similar reaction rates compared to pure CO 
ices. The presence of C0 2 in interstellar ices with CO there- 
fore does not affect the formation of H 2 CO and CH3OH. 

- Hydrogenation of CH3CHO leads for -20% to C 2 H 5 OH, 
showing for the first time that a thermal hydrogenation re- 
action can be responsible for the C 2 H50H abundances de- 
tected in dense interstellar clouds. Other reaction products 
are H 2 CO, CH3OH (15-50%) and CH 4 (-10%). Due to the 
energy released a fraction of the produced species may evap- 
orate into the gas phase upon formation. 
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